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Q ' In an investigation of the cause of the cosmic microwave background 

|H , decrement in the field of the z = 3.8 quasar pair PC1643+4631, we 

^ ' have carried out a study to photometricaUy estimate the redshifts of 

galaxies in deep multi-colour optical images of the field taken with the 
WHT. 

^> , To examine the possibility that a massive cluster of galaxies lies in 

^ • the field, we have attempted to recover simulated galaxies with intrinsic 

5^ , colours matching those of the model galaxies used in the photometric 

redshift estimation. We find that when such model galaxies are added 
to our images, there is considerable scatter of the recovered galaxy 
redshifts away from the model value; this scatter is larger than that 
expected from photometric errors and is the result of confusion, simply 
due to ground-based seeing, between objects in the field. 

We have also compared the likely efficiency of the photometric red- 
shift technique against the colour criteria used to select z > 3 galaxies 
via the strong colour signature of the Lyman-limit break. We find that 
these techniques may significantly underestimate the true surface den- 
sity of z > 3, due to confusion between the high-redshift galaxies and 
other objects near the line of sight. We argue that the actual surface 
density of z « 3 galaxies may be as much as 6 times greater than that 
estimated by previous ground-based studies, and note that this con- 
clusion is consistent with the surface density of high-redshift objects 
found in the HDF. Finally, we conclude that all ground-based deep 
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field surveys are inevitably affected by confusion, and note that reduc- 
ing the effective seeing in ground-based images will be of paramount 
importance in observing the distant universe. 



1 Introduction 

The deep five-colour photometry (Paper I) of PC1643-I-4631 allow an at- 
tempt to estimate the redshifts of objects in the field. To make the best use 
of our images, we must have an equivalent set of colours of galaxies derived 
from either observations or models. Composite galaxy spectra, such as those 
compiled by [ Coleman et al., 19801] (hereafter known as CWW), can be con- 



volved with the telescope transmission functions through each filter so that 
the colours of such a galaxy at various redshifts can be modelled, providing 
us with a "no-evolution mode" set of colours. The original CWW spectra 
cover wavelengths from 1500A to lOOOOA, while the filters effectively sample 
from 3000 A to 9000 A, so such an approach would provide model colours for 
z < 1. Extending these spectra further into the UV allows models for higher 
redshifts to be calculated and this can be achieved by careful matching of 



UV spectra, eg from the [Kinney et al., 1993] atlas, with the blue end of the 



optical spectra. The different morphological types of galaxies must also be 
taken into account - spiral galaxies tend, for example, to be considerably 
bluer than ellipticals - so different sample spectra must be compiled for each 
morphological group. 

The main restriction encountered with these composite spectra is the 
difficulty in adapting these spectra to account for evolution. This is not 
really a problem at low redshifts, but to extend these models back to the 
possible formation redshifts at zj = 5 one must include the evolution of these 
galaxies. The best one can do at present is to use the existing models, in 
which typically there is a burst of star formation at some formation redshift, 
Zf, followed by exponentially decreasing star- formation, and calculate the 
consequent colour evolution. 



2 Modelling the galaxy colours over < z < 4 

The composite spectra from actual galaxies provide the most represen- 
tative colours for the low redshift galaxies in the field (see, for example 



[ Gwyn, 199^ ]). At higher redshifts there are significant differences between 



the colours derived from empirical and simulated (in this case using the 
Bruzual &: Chariot (B&C hereafter, j Bruzual and Chariot, 199^ ) code) spec- 



tra: this is the result of including evolution in the simulated spectra. In 
order to achieve a consistent set of colours for each morphological type for 
redshifts between and 4, one cannot simply rely on either the B&C sim- 
ulations or the CWW spectra - the former do not match the low-redshift 
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galaxy spectra so well, while the latter do not include any evolution and can 
not, in any case, be used beyond z 1.5 since the UV end of the spectrum 
starts to move out of the U filter at this redshift. 

Colours derived from both sets of galaxy spectra are shown in Figure ||. 
Note that all the CWW spectra are much redder in both U — G and G — R 
at z = 1.5 than the B&C spectra. At this redshift, the star formation rate in 
the B&C spectra is producing considerably more flux in the UV than is seen 
in local galaxies, consistent with the apparent star formation history (see, 
eg. [|Madau et al., 199^ ). The result is that across the optical observing 
band (3000A - 8500A), the spectra of star-forming galaxies are close to 
flat spectrum for 1.2 < z < 2.3. The success of the Lyman-limit imaging 
searches for z > 3 galaxies using B&C-derived spectra indicates that at 
least some galaxies match these simulated spectra at high redshift. 

The appr oach adopted by Gwyn (|Gwyn, 1998 1) has changed since his 



1996 paper (| Gwyn and Hartwick, 1996|| ) on photometric redshifts in the 



Hubble Deep Field (HDF). Previously, he relied on using spectra derived 
entirely from B&C simulations. Now, where HDF galaxies are found to be 
at z < 1.5 using the B&C spectra, he derives new photometric redshifts 
from colours based on the CWW spectra and reports improved accuracy 
when compared against spectroscopic redshifts. However, this results in a 
large hiatus between the colours of a 2; = 1.6, B&C-derived galaxy and a 
z = 1.5, CWW-derived galaxy, which is unsatisfactory, as can be seen in 
Figure |l|. 

To attempt to make use of the CWW spectra to more accurately estimate 
the colours of low-redshift galaxies, we have used the following approach: 

for z < 0.5, galaxy colours are derived only from the 
CWW empirical spectra; 
for 0.5 < z < 1.5, galaxy colours are interpolated smoothly between the 
two sets of models, with colours at 2; ~ 0.6 being 
closest to the CWW galaxy colours, and the colours at 
z ~ 1.5 being closest to the B&C colours; 

for z > 1.5, galaxy colours are derived only from the 
B&C simulated spectra. 

The final galaxy models used are illustrated in Figure The CWW 
spectra are over-plotted to show the differences between the models. 

Four morphological types were used, corresponding to E/SO, Sbc, Scd 
and Irr, and spectra were produced at redshifts from 0.0 to 4.0 at 0.1 in- 
tervals. To cover the colour space more completely, these four models were 
interpolated linearly between the model types and in redshift, to create 13 
models at redshift intervals of 0.025. 
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3 Matching the model spectra against the obser- 
vations 

Since the number of counts received even in a small faint aperture is large due 
to the plates being background limited, one can treat all errors as gaussian 
about the mean flux. Fluxes (unnormalised) for the model galaxies are 
derived from the model spectra by convolution with the telescope response 
for each filter, including atmospheric absorption and the effect of scattering 
within the telescope light path. 

To compare the models against the observed fluxes, we used a test, 
with 



where is the observed flux in filter i, (Ji is the error in observed flux 
in filter i, and rrii is the model flux in fllter i. To obtain the optimum flt 
between the model and the observed fluxes, we minimised with respect 
to a. 

However, given the range of models available, it is important to deter- 
mine the range of best-fit models which are consistent with the observations, 
and for this we used a relative likelihood method. The likelihood of obtain- 
ing an observed set of fluxes fi with errors ui from model fluxes mi (where 
rrii is a function of redshift, z, and morphological type, T) is 



This assumes that the errors in the observed fluxes are described by a 
Gaussian distribution - a good assumption. From this likelihood distribu- 
tion, the modal redshift can be determined by identifying the model with 
the highest likelihood. 

Since one can apportion a probability to each model, one can use this to 
determine not only the most likely model, but also the mean redshift and 
the range of acceptable redshifts which are consistent with the observations 
within some relative confidence level. This is important, since for any set 
of models the distribution of likely models may not be mono-modal. We 
can use the difference between the mean and modal redshifts as a simple 
discriminator to identify cases where the redshift is not well constrained. 




(1) 





4 Redshift distribution 
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4.1 Overview 



The two redshift distribution histograms, based on mean and most probable 
redshift, are plotted in Figure |3|. The main features of both these diagrams 
are the two peaks: one at < 2; < 0.6, and the other at z ~ 2.4, with only a 
small number of objects lying in the range 1.0 < z < 2.2. This is reinforced 
by studying the distributions of objects in space against redshift ~ these are 
plotted in Figure |. 

It is clear that the apparent lack of objects in this redshift range is not a 
real effect - there cannot really be so few objects in this region. Identifying 
the cause of this effect is important for understanding the limitations of the 
photometric redshift technique. 

From the model colour tracks plotted against the actual objects colours 
(Figure it is clear that the close-packing of the model colours in the 
redshift range 1.0 < z < 2.1 is a problem. This arises because the continuum 
of the galaxy spectra is almost flat for all spectral types at wavelengths 
between 1400A and 2500A. Only the R — I colours provide any redshift 
sensitivity for objects with 1.0 < z < 1.5. The model colours effectively act 
as an 'attractor' for objects with similar colours - 'attracted' objects are 
labelled with the redshift and type of the model they end up with. Thus 
closely spaced models are less effective at distinguishing objects, unless there 
are significantly more objects concentrated in that area. Even then, the 
similarity of the objects is sufficient to make the accurate determination of 
the redshift impossible in this redshift range given these filters. The result 
is that objects with redshifts 1.0 < z < 2.1 end up being pushed to lower 
{z ~ 0.8) or higher (z ~ 2.4) estimated redshifts. 

Beyond the limitations imposed by the filter selection, it is also notable 
that the number of objects that are close to the models SED colours is small 
in the redshift range 1.0 < z < 2.0. This suggests that the model colours 
being applied here do not directly correspond to observed colours in the 
field. This raises several possibilities: 

(1) the model colours are inaccurate; 

(2) there are significant systematic errors in the catalogue calibration; 

(3) the model colours are correct, but the observed colours differ due to 
reddening caused by the presence of dust in these objects. 

The accuracy of a catalogue's photometric measurements can be qual- 
ified by examining the accuracy of the photometric calibration, both as 
determined from the original calibrators and in comparison with other pub- 
lished results. All the calibrators used in the creation of the catalogue have 
high signal-to-noise, and even including the possibility of poor background 
level determination, the calibrations are at worst inaccurate by ±0.05 mag- 
nitudes in all filters. This is not sufficient to bias the results as dramatically 
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as observed here. In comparison with other pubUshed photometric mea- 
surements for the field, such as Hu & Ridgway's multi-filter imaging around 
Quasar A ( [ [Hu and Ridgway, 1994| ] ) , the results are all consistent within the 
photometric errors. 

Since the model colours here are based on unreddened spectra, it seems 
reasonable to suggest that the lack of objects with such blue colours is due 
to reddening of these objects by dust. This would have the greatest effect in 
the U &i G colours, since these are the ones with the shortest wavelengths. 
If we examine the spread of colours of the objects m G — R and R — I 
(Figure ^) we can see there is better coverage of the model tracks m. R — I 
than va. G — R. The objects at R — I < —0.5 are all within 2-a of i? — / = 0; 
ie this is most probably photometric scatter. This is particularly noticeable 
in R — I because of the lower depth of the / image compared with the rest 
of the broadband colours. 



4.2 Most likely or mean redshift? 

Comparing the two estimates shows several interesting features, as illus- 
trated in Figure 0. In 1486 out of 1665 objects with R < 26.0 , there is less 
than 0.1 difference between the modal and mean redshift. 63 objects differ 
by 0.1 - 0.2 in redshift, and 116 differ by more than 0.2. For those objects 
whose redshift determinations differ by greater than 0.2, these divide into 
three major groups: those which have a modal redshift < modal z < 0.6 
but have a large mean redshift, a set with 2.5 < modal z < 3.0 but a mean 
redshift of 1.8 < mean z < 2.7, and a set with modal z ~ 4.0 but a lower 
mean redshift. 

In Figure ^, we plot the precision of the redshift estimates as a function of 
R magnitude. The objects which show the greatest differences are all faint 
{R > 24.0); the fitting is looser due to the greater flux uncertainties. 
This in turn leads to more models sharing similar model-fitting likelihoods 
and a greater spread of redshifts with which the models are statistically 
acceptable. Examination of the precision of the redshift estimates against 
the modal redshift (Figure ^ shows that the majority of objects with < 
modal z < 0.8 have small standard deviations, while the fitting is generally 
less constrained at higher redshifts {z > 2.0) although the average standard 
deviations even at z ~ 3 is still Az ~ 0.1. 

The low-redshift group (0 < modal z < 0.6) consists of 102 objects, 
all of which have 24 < R < 26. Some examples of the object spectra and 
redshift likelihood distributions are shown in Figure O, and an explanation 



of the information on these graphs precedes in Figure IC 

There appear to be two groups of objects which show significant differ- 
ences between mean and modal redshift: faint blue galaxies which, by the 
nature of their flat and relatively featureless continuum, are difficult to ac- 
curately determine the redshift for without good photometric constraints; 
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and objects with weak constraints on the U magnitude which can mimic 
both z « galaxies and spirals at z 1.0. 

The intermediate (2.5 < modal z < 3.0) and high (3.8 < modal z < 4.0) 
redshift galaxies also show a fairly flat spectrum in G,V,R & /, and red 



U — G > colours (see Figures O & 13 for examples). Again this ambiguity 



between a low redshift E/SO/Sa galaxy with a strong 4000A break and a 
high-redshift starforming galaxy with absorption due to the Ly— a forest 
and extinction below 912A causes problems for the photometric redshift 
models. J, H and other IR imaging would break this degeneracy, as would 
significantly deeper U imaging. 

In summary, the modal redshift appears to give the better redshift de- 
termination. When ambiguities arise in the determination of the redshift, 
the standard deviation of the redshift is sufficiently sensitive to mark the 
problem cases. 

4.3 Low— redshift galaxies — z < 1.0 

The accuracy of the photometric redshifts is likely to be best in this sub- 
sample, especially in the range < z < 0.5. In a blind test carried out by 
Hogg et al, the rms deviation of the photometric redshift from the spectro- 
scopic determination was routinely Az < 0.1, using similar methods and a 
similar range of wavelengths covered by the filters used here. 

Examples of redshift estimations are shown in Figures 0-0. These 
also include estimations of galaxy morphological type, made possible by 
the wide spread of colour space covered by the photometric models at low 
redshift. 

A few spectroscopic redshifts for some of the galaxies in this field are 
already available, obtained by one of us with the WHT allowing one to 
check the accuracy of the photometric predictions. Table |l] shows that 
the photometric redshifts are accurate to Az ~ 0.15, although the sys- 
tematic offset suggests that the models can be further improved. A sim- 
ilar overestimate is seen in [ Lanzetta et al., 1996(1 where the photometric 



models are based on redshifted CWW spectra alone, where the rms scatter 
was reported as Az ~ 0.15. Comparing the published redshift for HRIO 
( IGraham and Dey, 1996f| ) , it is clear that the reddening of the spectrum in 
this z > 1 galaxy has resulted in an erroneous estimate. 

4.4 Intermediate— redshift galaxies — 1.0 < z < 2.5 



As discussed in section ^jj there are too few objects identified in the range 
1.0 < z < 2.5. It is important to determine what redshifts objects in this 
redshift range have in fact been given. Assuming that the models are at least 
a reasonable pointer to the correct colours, possibly affected by reddening, 
objects at these redshifts will have been mis-classified towards those models 
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which have colours similar to the 1.0 < z < 2.1 models. The most likely 
cases are therefore the low-redshift Irregular galaxies, which have relatively 
flat spectra in the visible wavelengths, and those models at redshifts z ~ 0.8 
and z ~ 2.3 which are adjacent region we are interested in. This can actually 
be seen in the histogram of objects against modal redshifts (Figure ^) - there 
is an apparent excess of low-redshift objects at z = along with another 
peak in the distribution around z = 2.4. No excess is seen at z ^ 0.8 but 
any increase here would be likely to be a small fraction of the real objects 
at this redshift. 

Those objects which are identified atl.0<2;<2.5- examples are 



shown in Figure |18| - are only weakly distinguished in morphological type. 
Additionally there is a degeneracy between morphological type and redshift 
estimate for galaxies with redshifts 1.0 < z < 1.6 - at higher redshifts 
there are few major differences between the galaxy colour tracks and hence 
no significant discrimination between morphology based on colour evidence. 
To remove these degeneracies would require more IR imaging to detect the 
4000A break - this would also improve the accuracy of the redshift detections 
in this redshift range. 

4.5 High— redshift galaxies — 2.5 < z < 4.0 

In Paper I, we have tried to identify galaxies at z ~ 3 via the 912A break. 
Photometric redshifts should prove effective at determining which galaxies 
lie at these redshifts since they make use of all the photometric measure- 
ments together in all filters, rather than just the U,G & R. 



Examples are shown in Figure 19. It is not possible to differentiate 



between different morphological types at these redshifts - the colours are too 
similar between the model galaxies. This is partly a result of similar star- 
formation histories between the model galaxies at these redshifts; this would 
be invalid if, for example, there were several phases of galaxy formation. 

By selecting galaxies with 2.8 < modal z < 3.5, R < 25.5 and G > 2a, 
we get 46 Ly-break candidates, almost twice as many as with the colour- 
criteria. The distribution of these galaxies (see Figure ^) also (cf Paper I) 
shows an apparent hole at a point mid-way between the quasars, despite 
the increase in surface density of candidates. 



5 Clustering of objects 

We can make use of the redshift determinations to look at the correlation 
functions for redshift subsets of the catalogue. The results of these are 
plotted in Figure 

The lowest redshift range {0 < z < 0.5) has 558 objects, and there 
are 352 objects in the range 0.5 < z < 1.0. As a result, the statistics for 
these redshift ranges are sufficiently good to show that there is measurable 
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clustering in this field at low redshifts. Fitting the standard power-law 
curves (^^^6"^) to angular scales < ^ < 150 arcsec gives 0.82 and 

/? « 1.05 for 0.0 < z < 0.5, and « 1.88 and (3 ^ 1.74 for 0.5 < z < 1.0. 
These are plotted in Figures 22 & 23. Also plotted in these figures are the 
correlation functions for subsets of the sample ordered in redshift. While 
the noise in these subsets is worse and no strong conclusions can be drawn, 
the apparent clustering of the 0.7 < z < 0.9 objects at small angular scales 
is interesting. 

In the rest of the z bins there is little significant signal, with the least 
insignificant being 1.5-a clustering at 240" in 2.0 < z < 2.5. Necessarily, 
there is no coverage of 1.0 < z < 2.0. 



6 Artificial clusters 

We now test the ability of this sort of analysis to recover real features from 
the field. We use these redshift colour models to place a realistic spatial 
distribution of artificial galaxies into the images and examine the likelihood 
of detecting a cluster towards PC1643+4631 A&B using colour methods. 
We can find no similar experiment in the published literature on any field. 



6.1 Theoretical magnitude and spatial distributions 

To mimic a cluster of galaxies, one needs distributions for both the mag- 
nitudes and positions of the galaxies on the sky. For the former we as- 
summed a Schechter luminosity function ( |Schechter, 197(^ ]), while for the 
latter we assumed a Hubble distribution, where the surface density is de- 
scribed by c7o/(l + {R/Rc?) ( jWu and Hammer, 1993|| ). We used a core 
radius, Rc = 200kpc for the spatial distribution (Jones et al., 1997, found, 
for PC1643, a core radius of ^ 300kpc), and used half-light radii of lOkpc 
for the galaxies. 

To determine the magnitudes of galaxies in a cluster at some red- 
shift, we used the results from the Canada France Redshift Survey (CFRS) 
( jLilly et al., 1995[| ), because it is one of the best measured in a constant 
rest frame waveband and it provides us with the best-fit S-band Schechter 
luminosity function for 2; ^ 0.5 galaxies ([ [Lilly et al., 199"5[ ). The CFRS 
contains many objects at redshifts out to z ~ 1.1, and provides the best es- 
timate of the luminosity at 2: ~ 0.5. Additionally, for z < 0.9, the i?-band is 
redshifted to wavelengths covered by the U,G,V,R and I filters. Taking the 
values a = 1.03, absolute magnitude M^^{B) = —21.04, we determine the 
observed total magnitudes from these galaxies using (see jLilly et al., 1995f| ) 



Jab = ^^AB + 51og(L>L/10pc) - 2.51og(l + z) - {B - h), (3) 
Dl'is the luminosity distance, Iz is the I band shifted to shorter wavelengths 
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by a factor of (l + z), and B is the magnitude in the rest-frame B band. The 
value of {B — Iz) is taken from the photometric SED models, and hence the 
value of B must be calculated by interpolation between the C/, G and other 
broadband filters. From this, observed magnitudes can be calculated for the 
redshift range < 2 < 1.5. For z > 1.5, B moves to longer wavelengths 
than /. A simple extrapolation of the R — I colour can be used to obtain a 
first-order estimate of the rest frame B band magnitude. This should prove 
to be a reasonable approximation of the real value given that the 4000A 
break is weak in galaxies which are strongly star-forming, which is the case 
for the SED models used at z > 1. 

We do not attempt to model the evolution of the luminosity function. 
While this is a shortcoming of this method, at z < 0.9, the maximum likely 
error is ±0.5 mag (see [ Lilly et al., 1995[] ). Beyond 2; ~ 1, the evolution of 



the luminosity function is poorly constrained, but it is reassuring that the 
luminosity function of z ~ 3 galaxies have a luminosity function best fitted 
by -21 ( jPickinson, 1997[| ). 

At low redshifts, there are significant colour differences between differ- 
ent morphological galaxy types, with spirals having much bluer colours than 
ellipticals. While the ratio of ellipticals to spiral galaxies in general is about 
60:40, in the cluster environment the proportion of elliptical galaxies is of- 
ten much lower, closer to 40:60 (E+SO:S+Irr) ([ pressler et al., 1997| ]). To 
simulate this spread of colours, we opted to create frames with just one 
morphological type in each. To keep the processing as simple as possible 
two plates were made, with one having E/SO type galaxies and the other 
Scd/Irr types. 

6.2 Method 

The processes required to simulate each cluster at each redshift and make 
photometric measurements of the simulated cluster galaxies after superim- 



position onto the PC1643 are illustrated in Figure |2^. Simulated clusters 
were made at redshifts from 0.1 to 4.0 in steps of 0.1. 

We used the ARTGAL package in IRAF to create artificial clusters of 
galaxies with the required properties. For each 'cluster' only two images 
were made, one each for elliptical and spiral galaxies, made to have the 
same zero point and exposure time as the / image, and zero background 
level. No additional noise was simulated because the real images are all 
background limited. 

Recent publications discussing the distributions of the different morpho- 
logical types suggest that there is some segregation, with steeper velocity 



dispersions for later-type galaxies (see, eg [jAdami et al., 1998| ). We have 
not attempted to mimic these results here, since our aim is to examine 
the ability of the tests previously carried out on this field to recover cluster 
members, rather than to produce a completely accurate cluster distribution. 
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Having generated the I-band spiral image at each redshift, we used 
the Sbc model SEDs to determine the colours of the simulated spirals in 
R,V,G &: U, and IMARITH was used to scale the / image to mimic the 
correct colours in the other broadband images. We used a similar procedure 
to produce R,V,G & U images from the /-band elliptical image using the 
E/SO model SEDs. These I , R, V , G k, U spiral and elliptical images were 
summed together to create composite images, and a catalogue was made 
from the R composite image using FOCAS, which we will refer to as the 
matching composite catalogue (MCC). 

These composite images were then added onto the real PC1643 images. 
The resulting R image was processed using FOCAS to produce a catalogue 
based on R isophotal apertures: catalogues were then made m. I,V , G hU 
using the R isophotes and the five resulting catalogues were then matched 
using MATCH against the MCC. To extract the simulated galaxies from 
these newly created FOCAS catalogues, 'stripped' catalogues were made 
including objects only if they corresponded to the positions indicated by 
the MCC - ie these catalogues should contain only the information on the 
measureable simulated cluster galaxy members. These stripped catalogues 
were then given the same photometric redshift analysis as the real objects 
in the PC1643 field. 



6.3 Results 

Full-colour images showing the artificial cluster at various redshifts are 



shown in Figures 25 & 26, both in isolation and superimposed on the PC1643 
field. While it is easy to recognise the cluster on its own, it becomes im- 
possible to differentiate the cluster against the other objects in the field at 
z > 1.0, and the contrast is poor even at z = 0.6. 

More importantly, it is the ability of the catalogue creation and analysis 
software to detect a cluster which is of greatest interest here, particularly at 
faint magnitudes R > 24. We have plotted the U — G vs G — R colours for 



each of the simulation images illustrated in Figures ^ Even at z = 0.2 
the fainter members of the simulated galaxies show significant photometric 
scatter away from the expected colours, and this is reflected in the estimates 
of modal redshift: there is a peak at z = 0.2 as expected, but over a third of 
the sample are more than 0.1 away from this redshift, and 9 of the simulated 
galaxies are mistaken for high-redshift galaxies at z > 3. 

At z < 1.0, the photometric-redshift technique still produces redshift 
estimates around the expected value, with increased scatter as the simulated 



cluster becomes fainter with increasing redshift (see Figures 29 & ^0|). This is 
also reflected in the increased scatter of the U — G and G — R colours, which, 
it should be stressed, is often greater than S-o", where a is the photometric 
error based on the noise for the object. This is the result of pollution of 
the isophotal apertures by nearby neighbours, which is inevitable in a deep. 
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ground-based field like this one. Given that the density of objects in the 
field including all detections down to the surface brightness limits of 3-0" in 
an area of six pixels is approximately one object in every 5" x 5" box, and 
that the average isophotal area of a object with R = 25.0 is some 5 arcmin^, 
this means that at least one in five faint objects will overlap another faint 
object, so that in this set of deep images, objects are being confused. 

To determine whether a significantly brighter cluster of galaxies would 
be visible in this field, we repeated the simulation with all the magnitudes 
brightened by one magnitude. Bearing in mind that this is the most extreme 
error likely in the simulation, this is also a good test of the sensitivity of 
these simulations to the real magnitude distribution. Even with many more 
members, such a cluster would still be difficult to distinguish using colour- 
criteria or by excess surface density. 

At intermediate redshifts, 1.0 < z < 2.0, the photometric redshift tech- 
nique fails to correctly determine the redshift in almost all cases. Despite 
starting from simulated galaxies with the correct colours, the colours of the 
galaxies as recovered from the field are systematically scattered to redder 
colours. There may be a systematic effect resulting from selecting objects 
based on their R magnitudes — for example, objects which are polluted to 
fainter magnitudes in R are missed from the sample, whereas objects which 
recover fainter G magnitudes than originally simulated appear redder than 
they really are. There may also be a systematic effect from the different sur- 
face densities of objects in R and G at similar magnitudes - there are 763 
objects with 25 < -R < 26, and only 537 with 25 < G < 26 - so objects with 
G — ~ stand proportionally more chance of being brightened in R than 
in G. Additionally, if the position of a simulated galaxy does coincide with 
that of a real object, it is likely that the real object has G — R > 0.5, and 
hence the pollution of the simulated object will result in redder measured 
colours. 

This failure to accurately detect galaxies in the range 1.0 < z < 2.0 
is a function of the filters used. Because there are no strong continuum 
features over the wavelengths spanned by the optical filters used here at 
these redshifts (approximately 1200A - 4000A) the photometric redshift 
are poorly determined. Additionally, because of the scatter of the objects 
colours away from the simulated colours, objects with redshifts in the range 
1.0 < z < 2.0 are scattered to lower and higher redshifts. This is clearly 
seen in Figures 29 & The result is that any histogram of galaxies over 
the redshifts estimated from photometric techniques using a similar range of 
filters to those used here will show a deficit of galaxies with 1.0 < z < 2.0 and 
peaks above and below this range against the true population of galaxies. 
Precisely this behaviour is seen in the redshift distribution calculated from 
the HDF four-colour photometry (in F300W, F450W, F606W and F814W) 
using photometric redshifts (see |Gwyn and Hartwick, 1996 1 and to a lesser 
extent [ gawicki et al., 1997 ]). The deficit is not quite so marked as it is in 
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the photometric redshifts estimates of the field of PC1643, due to sHght 
differences in modelling but the effect is still the same. It is notable that 
the photometric redshift estimates of [ Lanzetta et al., 1996| ] (LYFS) do not 
suffer from this effect, and appear to produce a smooth redshift histogram: 
the photometric models used by LYFS are based entirely on the CWW 
spectra extended into the UV without including the effects of evolution, 
with the effects of the Ly-a forest and 912A break taken into account. 
However, the galaxies at z f« 3 are predicted to have C/300 — V450 ~ 5 and 
V450 — -R606 1.0 by LYFS, which is considerably redder in [/300 than 
those detected by [ Steidel et al., 1996a ] which have U300 — V450 ~ 2 and 
V450 — -Reoe ~ 0.5, and therefore casts doubt onto their photometric models. 

To test just how much brighter a distant cluster would have to be to be 
clearly evident, we increased the magnitudes of every member of the simu- 
lated cluster ai z = 2.0 until the cluster just became evident in the images. 
We emphasize that significantly more than two magnitudes of brightening 
would be required to make the cluster stand out from the field. We stressed 
that this is far more than the largest expected error in the brightness of 
these simulations. 

At z > 2.3, the 912A break starts to extinguish a [/-band observa- 
tion, and by z ~ 3 provides a clear colour signature, ie ?7 — G > 2.0. 
This feature has been successfully used to detect high-redshift galaxies 
both with custom filters and in the HDF (eg [ Steidel and Hamilton, 1993 ], 
[ Steidel et al., 1996a|). Most of the galaxies previously detected at z « 3 
have R > 24.5 [ Steidel et al., 1996b [, which is consistent with the brightest 
members of the simulated cluster at the same redshift (Figure 31). Intrigu- 
ingly, while the photometric redshift estimations do peak at around their 
expected value, there is also a significant proportion which are mistaken for 
low-redshift galaxies (again, see Figures |2| & pOl). From the 123 simulated 
galaxies with colours consistent with 2.5 < z < 3.0 and measured R < 26.0, 
52 are identified as having z > 2.5 and 70 are identified as low redshift 
galaxies with z < 0.8. If this mimics the real situation for detecting 2; ~ 3.0 
galaxies, then it suggests that a significant number of high-redshift galaxies 
are missed if redshifts are estimated using photometric techniques. Note 
also that the limiting magnitudes of [ Steidel and Hamilton, 1993| ] etc. are 
similar to those of the PC1643 images. 

To investigate the efficiency of selecting z w 3 galaxies using the colour- 
criteria described in Paper I, we concatenated the simulated catalogues with 
2.5 < z < 3.5 and selected objects with R < 25.5 and G > 2a. While the 
simulations were carried out to mimic a cluster of galaxies, there is negligible 
overlap between the individual simulated galaxies in each catalogue. There- 
fore, while there is clustering of the galaxies in the field, it is important to 
note that the photometric spread and the results are similar to that which 
would be obtained for a random uniform spatial distribution. 

Figure 32 shows the measured U—G vs G—R colours of these 104 objects. 
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We note that the simulated objects are scattered from the simulated colours 
to such an extent that only 12 of the objects show measured U — G and 
G — R colours consistent with objects dX z ?>. Of these same 104 objects, 
48 are identified as being at 2 > 2 by the photometric redshift estimator 
(see Figure ^). Even if we take just the catalogue for the simulated cluster 
at z = 3.0, of the 12 objects with R < 25.5 and G > 2a, only two meet 
the colour-criteria, whereas eight are identified as having estimated z > 2. 
Therefore, five-colour photometric redshifts appear to provide a more effec- 
tive method of detecting high-redshift galaxies than the three-colour U GR 
method used by Steidel et al. This is not surprising since the extra photo- 
metric measurements provide extra constraints, which reduce the effects of 
individual photometric measurement errors. 

Comparing the angular sizes of the simulated objects against Ly-break 
galaxies observed in the Hubble Deep Field suggests that these galaxies 
should be about 2-3kpc in size ( |Dickinson, 199'^ ]), which is considerably 
smaller than the sizes used in this simulation, which are « lOkpc in size. We 
repeated the z = 3.0 simulation to assess whether the angular sizes used had 
a significant effect, with the angular sizes of the simulated galaxies reduced 
by a factor of four. This small angular size simulation had 15 objects with 
R < 25.5 and G > 2cr, of which three fulfilled the colour-criteria, whereas 
using photometric redshifts, 8 were identified at z > 2. Figure ^ compares 
the U — G vs G — R colours for the original and revised simulation. 

The increase in numbers of galaxies detected is due to the increase in 
the peak flux of galaxies due to smaller half-light radii - this makes it 
easier for FOCAS to detect these galaxies as their peak fluxes are more 
likely to be greater than the 3a detection threshold. It is important to note 
that decreasing the angular size of the simulated galaxies has not had a 
significant effect on the fraction of galaxies which fulfill the colour criteria 
for z ~ 3 galaxies: one might expect reducing the angular size to reduce 
the level of confusion observed in these simulations; however, while the peak 
flux increases and the area covered by the object is decreased for all objects, 
this also promotes simulated objects which were previously too faint to be 
detected above the detection threshold and it is these faint objects which 
are most likely to be significantly affected by confusion with the many real 
faint objects in the field. The effect of confusion on objects with faint (ie 
R > 24) magnitudes is to "average" the simulated object's colours with 
that of the confusing object, with the result that a large proportion of the 
sample is moved towards the average colours of the real objects in the field, 
ie U - G ^ 1.2 and G - R^ 1. 

From these results, it appears that the colour-criteria of 
[^teidel et al., 1995 ] is missing a large fraction, 5 out of 6, of the real 
population of z ~ 3 galaxies. This is the result of confusion due to the 
~ 1" seeing in the images; one would therefore expect to find considerably 
greater numbers of Ly-break galaxies in the HDF. There are 26 galaxies in 
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the HDF which have spectroscopic redshifts z > 2 (| Dickinson, 199^ ]), of 
which 12 have 2.8 < z < 3.5 and Vgoe < 25.5 (this is consistent with the 
spectroscopic redshifts of Ly-break candidates; see [ Steidel et al., 1996b ] 
and the faint magnitude hmit). The area of the Hubble Deep Field is 
4.5 arcmin^ and one would expect, given the published surface-density of 
"robust" Ly-break candidates of « 0.7 arcmin^ -2 ( [^teidel et al., 1998|| ) to 
find only 3 Ly-break candidates in the HDF. Because the HDF covers such 
a small field of view, there may be significant differences in the numbers of 
Ly-break galaxies between adjacent fields, but it is striking that the factor 
of apparent overdensity of Ly-break galaxies seen in the HDF is equivalent 
to the factor of Ly-break galaxies lost to confusion in the ground-based 
images. 

In [ Dickinson, 1997 1, the luminosity function of z ~ 3 galaxies is illus- 
trated, but the HDF count of z ~ 3 galaxies have been renormalised to the 
numbers seen in the ground-based images. We suggest that the surface den- 
sity seen in the HDF is a better estimate of the real surface density of z « 3 
galaxies, and there are approximately 6 times more z ~ 3 galaxies than cur- 
rently thought. This has important consequences for the understanding of 
the star-formation history of the universe and would suggest that the peak 
rate of star-formation occurred at an earlier epoch than currently thought 
(see [ [Madau et al., 1998[ ). 

To definitively identify the majority of Ly-break candidates and de- 
termine the real surface-density of these galaxies, observations with much 
improved seeing will be necessary; either by further HST observations to 
increase the area of sky observed and thereby reduce the effects of cosmic 
variance, or by deep ground-based observations with adaptic optics to re- 
duce the point -spread-function in the images. The effects of confusion can 
be quantified using similar simulations to those demonstrated here. 



7 Conclusions 

In using the multicolour photometric redshift technique on the objects found 
in the deep imaging of PC1643-I-4631, we have shown the following. 

(1) The technique is efficient at identifying objects at z < 1 where the 
photometric errors are small, ie at i? < 24, because the fitting is well 
constrained because the colours change strongly with redshift. 

(2) These photometric redshifts are in broad agreement (Az ~ 0.15) with 
the 4 spectroscopic redshifts in the range 0.62 < z < 0.81 in the 
PC1643-I-4631 field, whereas HRIO at z = 1.44 has a colour-estimated 
redshift of z = 0.65, indicating that either significant reddening can dis- 
tort the redshift determination or that at 1 < 2; < 2 the photometric 
technique fails. 
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(3) Ground-based morphological classification is possible on low-redshift 
(-2 ^ 0-7) objects with high signal-to-noise photometry (ie R < 24 
in the deep PC1643 images), although this work does not investigate 
its accuracy, due to the difficulties in obtaining secure morphological 
classifications for a significant number of galaxies in this field by manual 
or automatic means. 

(4) There are degeneracies between the morphological type and redshift in 
objects with 1.0 < z < 1.6, due to the models with different morpholo- 
gies having similar colours at differing redshifts. 

(5) No constraints can be placed on the morphology of high redshift (z > 
3) galaxies using these models since there is no significant difference 
between the colours of the model galaxies with different morphologies 
at z > 2. 

(6) 46 Ly-break candidates have been identified in the PC1643 images with 
2.8 < z < 3.5, R < 25.5 and detected at > 2a in G; the distribution 
of these candidates still shows an underdensity mid-way between the 
quasars. 

(7) The histogram of photometric redshift estimates shows two peaks, at z ~ 
0.7 and z « 2.3, and few objects with 1.0 < z < 2.0. This is probably 
an artifact of the photometric redshift technique, arising from the close- 
spacing of the model colours at these redshifts, and the histogram is 
therefore not indicative of the real distribution of redshifts. 

(8) By calculating both modal and mean photometric redshifts, as well as 
a standard deviation for the likelihood distribution, for each object, 
ambiguous photometric redshift estimates can easily be flagged using 
|modal z — mean z\ > 0.2, giving an advantage over previous modal- 
only methods. 

(9) There is evidence of clustering in the field at low {z < 0.9) redshift which 
is consistent with other published results. 

We have carried out simulations to examine the ability of the photomet- 
ric redshift technique to accurately recover simulated galaxies at various 
redshifts. These show the following. 

(10) For simulated galaxies at 0.1 < z < 0.3, faint E/SO galaxies are mis- 
takenly identified as z > 3 galaxies because of the red U — G colours. 
An E/SO with R ^ 2A will often have no detected flux in ?7, and be 
mistaken for a Ly-break candidate. 

(11) For simulated galaxies at 0.1 < z < 1.0, the average estimate of red- 
shift is consistent with the simulated redshift, but there is a spread of 
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Az ~ 0.2 in the distribution. This is partly due to the photometric er- 
rors causing the colours of the object to move away from the simulated 
colours, but also due to contamination of the object colours by other 
objects in the field. 

(12) For simulated galaxies at 1.0 < z < 2, the photometric redshift is un- 
reliable and misclassifies these galaxies to lower (z ~ 0.6) and higher 
(z ~ 2.2) redshifts, resulting in peaks in the photometric redshift dis- 
tribution which are not representative of the real redshift distribu- 
tion. This is a fundamental problem with photometric redshift esti- 
mates based only on optical imaging. A similar problem is seen in 
several publications in the literature but has not been appreciated (see 
|Gwyn and Hartwick, 1996 1 and [ Sawicki et al., 1997| ], who both have a 



decrease in numbers of galaxies between 1 < z < 2 but considered it 
a real observation), which has resulted in erroneous conclusions being 
drawn. 

(13) For simulated galaxies at z > 2.5, the photometric technique is effective 
in detecting about half the sample as being at z > 2.5, with the rest of 
the sample being mistaken for low redshift galaxies. 

(14) Photometric estimates may prove effective in selecting high redshift 
galaxies at z ~ 3, where the U — G > 2 colours and fairly flat spec- 
trum in G,V,R and / provide additional constraints to the three colour 
U,G,R method currently used for successfully selecting z ^ 3 galaxies. 

(15) The scatter of the simulated galaxies away from their simulated colours 
is two to three times greater than that expected from the photometric 
errors alone and arises as a result of confusion between the simulated 
galaxies and the real objects in the field. Since the amount of confusion 
is directly related to the resolution of the images, the WHT images will 
be far more affected by confusion than equivalently deep HST images. 

(16) Because of the effects of confusion, the current published estimates of the 
surface densities of Ly-break galaxies based on ground-based imaging 
are under-estimates by a factor of ~ 6; such an increase in the surface 
density of these galaxies is consistent with the surface densities seen in 
the HDF. 

(17) Any third image of a single quasar causing the PC1643 pair is extremely 
likely to be confused and therefore unrecognisable. 

(18) Deep ground-based CCD imaging without adaptive optics is inevitably 
confusion-limited at R > 26.0. Imaging with smaller point -spread- 
functions is the only way to alleviate this. 
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Simulations have been carried out to investigate the difficulties in iden- 
tifying a cluster at various redshifts in the PC1643 field. These show the 
following. 

(19) Even with our many optical colours, it is difficult to visually identify a 
cluster of galaxies at z ^ 0.4 against the other objects in any deep field. 
A cluster at z = 1 with galaxy luminosities similar to those seen locally 
would be lost against the other objects in the field. 

(20) Even a very rich cluster at 2; 1.0 would be impossible to distinguish 
from the field on the basis of this U, G, V, R and I imaging and hence 

the evident absence of a cluster in these l.l" resolution images is entirely 
consistent with a distant cluster producing the CMB decrement seen. 

(21) A cluster at z = 1.0 in which each member was 1 magnitude brighter 
than expected would still be difficult to detect; however, if the members 
were more than 2 magnitudes brighter, the cluster would be evident in 
the images. 
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Spectroscopic Photometric 
Object redshift redshift 



1 
2 
3 

4 

HRIO 



0.659 ± 0.002 
0.670 ± 0.002 
0.670 ± 0.002 
0.810 ±0.003 
1.44 ±0.01 



0.55 ± 0.03 
0.57 ± 0.04 
0.52 ± 0.05 
0.62 ± 0.08 
0.65 ± 0.08 



Table 1: Comparison of spectroscopic and photometric redshifts in PC1643 
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Figure 1: Comparison of U-G and G-R colours based on simulated galaxy 
spectra for four morphological types against the empirical galaxy spectra 
taken from Coleman, Wu and Weedman. The points are at intervals of 0.1 
in redshift 
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Final model colours 
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Figure 3: Distributions based on photometric redshifts for all objects with 
R < 26.0 
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Figure 4: Angular distributions versus redshift for all objects with R < 26.0. 
EW position is the offset from the centre of the field along the X-axis of 
the CCD. NS position is the offset form the centre of the position along the 
Y-axis of the CCD. 
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G — R /magnitudes 

Figure 5: U-G vs G-R plot overlaid with the revised models colour tracks 
used in the determination of the photometric redshifts. 
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Figure 6: G — Rvs R — I plot overlaid with the revised models colour tracks 
used in the determination of the photometric redshifts. 
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Figure 7: Most likely redshift plotted against the mean redshift 
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Precision of redsliift estimates 

All objects have R<26.0 
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Figure 8: Standard deviations of the redshift distributions for each object 
plotted against R magnitude. 
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Precision of redsliift estimates 

All objects have R<26.0 
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Figure 9: Standard deviations of the redshift distributions for each object 
plotted against modal redshift. 
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Figure 10: Explanation of the layout of the likelihood distribution graphs 
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Figure 11: Examples of five objects with low modal redshifts but complex 
likelihood distributions 
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Figure 12: Two intermediate modal redshift objects with complex likelihood 
distributions 
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Figure 13: Two high modal redshift objects with complex likelihood distri- 
butions 
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Figure 14: Two low-redshift E/SO-type galaxies 
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Figure 15: Two low~redshift Sbc-type galaxies 
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Figure 16: Two low-redshift Scd-type galaxies 
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Figure 17: Two low-redshift Irr-type galaxies 
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Figure 18: Six examples of galaxies with fits to model spectra 1.0 < z < 2.5. 
Note that there is httle constraint on the morphology of the objects shown 
here. Additionally, at z ~ 1.3 there is a degeneracy between morphology 
and redshift. 
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Source no.; 2514 Mean redshift: 3.4 redshift s.d.: 0.23 
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Figure 19: Examples of galaxies with fits to model spectra z > 2.5. There is 
no possibility of morphological classification at these redshifts. Note that at 
z ^ 3.4 the effect of the 912A break moving into the blue end of the G band 
starts to redden the G — R colours, and that this feature is identified by the 
photometric redshift technique more strongly than the faint U magnitudes, 
which have large photometric uncertainties. 
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PC1 643+4631 - high redshift candidates 
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Figure 20: The distribution of the Ly-break candidates, which are marked 
as circles, selected using R < 25.5, G > 2a and 2.8 < modal z < 3.5. 
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2-point correlation functions 

Sub-samples based on modal redshift ranges 
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Figure 21: Comparison of the clustering of the galaxies in the field in differ- 
ent redshift ranges 
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Figure 22: Correlation function for galaxies with modal redshifts < z < 0.5 
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0.5<z<1.0 




44 



Simulated clusters of galaxies; 
Image and catalogue creation 
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Figure 24: Flow diagram showing the procedures used for simulating images 
of clusters of galaxies and obtaining a catalogue from those images 
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This figure is avaliable at 
ftp. mrao. cam. ac.uk:/pub/PC1643/paper2. figure50.eps 



Figure 25: Colour images of the simulated cluster at z = 0.4, z = 0.6 and 
z = 0.8. Images on the left are on a blank background; those on the right are 
superimposed on the PC1643 field. All these images use the same colouring 
scheme used in Cotter et al., 1998. 
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This figure is avaliable at 
ftp . mr ao . cam . ac . uk : /pub/PC 1 643 /pap er 2 . figure5 1 . eps 



Figure 26: Colour images of the simulated cluster at z = 1.0, z = 2.0 and 
z = 4.0. Images on the left are on a blank background; those on the right are 
superimposed on the PC1643 field. All these images use the same colouring 
scheme used in Cotter et al., 1998. 
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Figure 27: The histogram of modal redshifts as estimated from photometric 
measurements of the simulated galaxies. These graphs illustrate the effect 
of attempting to recover the simulated galaxies from the PC1643 field using 
FOCAS at z = 0.2, z = 0.4, z = 0.6 and z = 0.8. Only simulated galaxies 
with measured i? < 26.0 are shown. 
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Figure 28: The histogram of modal redshifts as estimated from photometric 
measurements of the simulated galaxies. These graphs illustrate the effect 
of attempting to recover the simulated galaxies from the PC1643 field using 
FOCAS at z = 1.0, z = 2.0, z = 3.0 and z = 4.0. Only simulated galaxies 
with measured i? < 26.0 are shown. 
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Simulated clusters of galaxies 

Recovery of redshift of a simulated cluster using photometric redshifts 
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Figure 29: Modal redshift estimates based on photometry from recovering 

the member galaxies of simulated clusters at z = 0.1 to z = 4.0 in steps of 
0.1 in z using FOCAS. Modal redshifts are shown for all simulated galaxies 
which have a measured i? < 26.0. 
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Simulated clusters of galaxies 

Recovery of redshift of a simulated cluster using photometric redshifts 
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Figure 30: Mean redshift estimates based on photometry from recovering 

the member galaxies of simulated clusters at z = 0.1 to z = 4.0 in steps of 
0.1 in z using FOCAS. Mean redshifts are shown for all simulated galaxies 
which have a measured i? < 26.0. 
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Galaxy counts in R 
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Figure 31: R magnitude distribution for the simulated cluster at z = 3.0 



52 




-10 12 



G — R /magnitudes 

Figure 32: U — G vs G — R colour-colour graph for simulated galaxies with 
2.5 < z < 3.5. All objects shown have measured R < 25.5 and G > 2-cr. 
The triangles denote l~a lower limits in U — G. The cross are stars taken 
from the Gunn &; Stryker database, and the dot-dash line is the bound for 
selecting high redshift candidates. The tracks plotted on this graph are the 
model colours for galaxies at various redshifts: the point at U — G = 2.2 and 
G — R = 0.3 is at z = 3.0, and points at larger U — G are at higher redshift 
in steps of Az = 0.1. 
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Galaxy counts in Modal Z 
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Figure 33: Histogram of the photometric redshift distribution measured from 
a simulated catalogue containing objects with simulated colours consistent 
with galaxies at 2.5 < z < 3.5. 
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Figure 34: Comparison of the U — G vs G — R colours between simulations 
with galaxies with angular sizes ~ lOkpc (left) and galaxies with angular 
sizes of ~ 2.5kpc (right). The dashed line denotes the limit of the colour- 
criteria used to select Ly-break galaxies. All objects plotted have R < 25.5 
and G > 2a: circles represent > la detections in U,G and R; triangles 
denote 1-a lower limits for U — G. 
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